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ABSTRACT 


The  goal  of  this  research  was  to  experimentally  demon- 
strate the  feasibility  of  a  cylindrical  shock  wave  to  serve 
as  the  major  pressure  recovery  mechanism  for  either 
cylindrical  gas  dynamic  lasers  or  chemical  lasers.   An 
analysis  was  conducted  to  determine  the  minimum  distance 
between  the  nozzle  exit  cylinder  and  the  shock  wave  location, 
r  -r  ,  required  for  power  extraction.   A  cylindrical  nozzle 
array  was  designed  using  the  axisymmetric  method  of  charac- 
teristics for  a  Mach  number  of  4.5  with  a  flexibility  of 
other  Mach  numbers.   A  chamber  incorporating  end  wall  boundary 
layer  bleed  and  ejectors  for  back  pressure  control  was 
designed  and  built.   The  apparatus  was  run  up  to  a  pressure 
of  18  atmospheres.   Three  flow  regimes  were  identified.   A 
value  of  (r -r  )/(r  -r  )  of  0.3  was  attained,  where  r  is 
the  shock  wave  location.   Flow  in  the  end  wall  nozzles 
consisted  of  a  Coanda  jet  with  a  subsonic  flow  region 
adjacent  to  the  wall;  flow  was  not  altered  by  boundary  layer 
bleed.   Except  for  the  end  wall  nozzles  a  cylindrical  shock 
wave  was  established. 
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I.   INTRODUCTION 

At  present  there  are  three  types  of  lasers  that  employ 
supersonic  flow:   chemical,  CO  electrical,  and  gas  dynamic. 
All  three  are  capable  of  producing  large  amounts  of  lasing 
power,  but  are  restricted  to  use  in  platforms  that  are  not 
limited  by  weight  and  volume.   Due  primarily  to  the  large 
diff users  necessary  to  restore  the  flow  to  ambient  pressure, 
large  volume  and  weight  are  required  to  handle  a  given  flow. 

A  cylindrical  geometry  gas  dynamic  laser  has  advantages 
not  found  in  the  other  supersonic  lasers  mentioned: 

(1)  diffuserless  operation 

(2)  elimination  of  nozzle  bending  stresses 

(3)  elimination  of  boundary  layer  corner  interactions. 

Such  a  geometry  would  result  in  a  compact  system  achieving 
lasing  by  rapid  expansion  in  the  radial  direction  rather 
than  linearly  as  is  presently  done.   The  supersonic  radial 
flow  is  then  slowed  to  subsonic  and  raised  to  near  ambient 
pressures  by  traveling  through  a  cylindrical  shock  wave 
sufficiently  downstream  of  the  nozzles  to  allow  for  a 
reasonably  large  lasing  cavity.   A  cylindrical  shock  wave 
eliminates  the  need  for  diffusers,  thereby  reducing  the 
weight  and  volume  of  the  device  considerably. 

Supersonic  conditions  are  achieved  from  flow  traveling 
through  axially  symmetric  nozzle  rings,  thus  eliminating 
nozzle  bending  stresses  found  in  linear  gas  dynamic  lasers. 
A  schematic  of  a  cylindrical  GDL  is  shown  in  figure  (1) . 
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Figure  (2)  shows  schematic  diagrams  of  two  types  of 
linear  GDL's.   Boundary  layer  buildup  on  the  four  walls  of 
the  lasing  cavity  interact  at  the  corners,  creating  increased 
possibilities  of  flow  separation  and  non-homogeneous  conditions 
in  the  lasing  medium.   By  having  a  flat  plate  at  each  end  of 
the  nozzle  stack,  the  cylindrical  geometry  would  eliminate 
such  interaction.   Since  the  radial  flow  experiences 
increasing  area  as  it  moves  outward,  the  boundary  layer  buildup 
on  each  plate  is  less  than  in  the  linear  case. 

A  cold  flow  test  of  a  number  of  nozzle  rings  was  conducted 
to  establish  a  cylindrical  shock  wave  and  to  determine  whether 
such  an  apparatus  is  feasible  for  use  in  a  cylindrical  gas 
dynamic  laser. 

A  theoretical  power  output  example  is  included  in  section 
II  of  this  report.   Section  III  contains  the  development  of 
the  nozzle  feed,  the  nozzle  ring  geometry  and  the  ejector  and 
boundary  layer  bleed  systems.   Experimental  results  are 
described  in  section  IV. 
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FIGURE  (2)   SCHEMATICS  OF  LINEAR  GDL ' S 
(FROM  REF.  5) 
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II.   CYLINDRICAL  LASER  GEOMETRY 

A.   PREVIOUS  GDL'S 

Linear  gas  dynamic  lasers  are  currently  in  use  and  have 
successfully  produced  laser  beams  in  the  order  of  a  megawatt 
of  power.   These  devices  are  quite  large  and  heavy,  thus 
limiting  their  use  in  operating  environments  (such  as  an 
aircraft)  where  such  parameters  are  limited.   These  limita- 
tions are  due  primarily  to  the  large  diffusers  necessary  to 
return  the  supersonic  flow  to  ambient  pressure. 

A  cylindrical  GDL  could  reduce  the  weight  and  volume 
necessary  to  handle  a  given  flow  by  expansion  to  supersonic 
velocities  in  the  radial  direction  and  returning  to  ambient 
pressure  by  means  of  a  cylindrical  shock  wave.   Other  advan- 
tages include  elimination  of  bending  stresses  and  reduced 
boundary  layer  interferences. 

Ortwerth  [Ref .  (1) ]  designed  a  cylindrical  source  flow 
device,  but  retained  diffusers  similar  to  those  used  in 
linear  GDL's.   The  diffusers  were  aligned  radially  as  shown 
in  figure  (3) .   Such  a  design  reduces  the  volume  required  of 
such  a  device  by  having  an  annular  lasing  cavity  rather  than 
a  linear  one.   Retaining  the  diffusers  does  not  allow  for 
the  substantial  reductions  in  weight  and  volume  that  a 
cylindrical  shock  wave  would  facilitate. 

Sponsored  by  the  Air  Force  Weapons  Laboratory,  the  AVCO 
Everett  Research  Laboratory  [Ref,  (2)]  conducted  an  extensive 
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theoretical  program  to  establish  the  feasibility  of 
cylindrical  source  flow  for  a  cylindrical  GDL.   Included 
in  the  program  was  a  cold  flow  test  of  one  nozzle  ring  which 
established  supersonic  flow  at  the  predicted  mach  number  but 
failed  to  create  a  cylindrical  shock  wave.   Results  showed 
that  homogeneity  of  the  lasing  medium  can  be  made 
acceptable.   Separation  of  the  boundary  layer  interfered  with 
the  establishment  of  a  cylindrical  shock  wave,  and  it  was 
recommended  that  some  type  of  boundary  layer  control  be  used. 

Figure  (4)  shows  the  assembly  used  in  the  AVCO  Everett 
test. 

B.   THEORETICAL  POWER  OUTPUT 

In  a  linear  GDL  the  diffusers  can  be  placed  far  enough 
downstream  so  that  a  large  lasing  cavity  is  assured.   In 
a  cylindrical  GDL  the  width  of  the  lasing  cavity  is  deter- 
mined by  the  shock  location.   To  extract  all  the  power 
available  from  a  given  flow,  the  cylindrical  shock  must  not 
occur  too  near  the  nozzle  exit.   Figure  (5)  shows  the  nozzle 
stack  with  the  lasing  mirrors. 

Assuming  a  Gaussian  intensity  distribution  on  each  mirror, 
the  intensity  I  is  given  by 


j<r>  =  we~r2°  •  (1) 
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Satisfying  the  boundary  conditions  and  re-defining  r  yields 


I(r)  =  I  e   re"r       kilowfts  (2) 

cm 


where  r  and  r  are  shown  in  figure  (A-l)  in  Appendix  A. 

Integrating  over  the  mirror  area  gives  the  total  power 


P  =  2  I     /    r  expt^-V^  ]   dr  (3) 

p  r  r  -r 

e  e 

An  alternate  method  of  finding  the  power  available  from 
the  flow  is 


Power  =  (mass  flow  rate)  x  (specific  energy) 


or 


P  =  2  re  Lpoao(p/po)  (ae/aQ)  Mq   x  (E)  kilowatts        (4) 


where  subscripts  "e"  and  "o"  refer  to  exit  and  stagnation 
conditions  respectively.   Equating  equations  (3)  and  (4)  gives 

rm         - 
lp   /   r  exp  [£l£-  V2~  ]2  dr  =  re  Lpoao(p/p0)  (ae/aQ)  Me  E 

e  e  (5) 


where  the  only  unknown  is  the  outer  radius  of  the  mirror, 
m 
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Performing  the  integral  in  equation  (5)  and  solving  for 

r_ 


m 


^£77 


rm  =Vr   +  c  (6) 

m   "  e 


4re  LPoao(a/ao)Me  E(P/P0> 
where   c  = 


^ 


£  Ip  erf  v/2) 


Using  the  values  for  the  constants  listed  in  Appendix  A, 
equation  (6)  reduces  to 


r   =V2.13  +  .2692  L  (7) 


where  L  is  the  overall  length  of  the  nozzle  stack.  Figure 
(6)  shows  the  minimum  shock  location  for  increasing  nozzle 
stack  length. 

For  six  nozzle  rings,  the  outer  radius  of  the  mirrors 
is  1.77  inches.   Using  this  value  in  equation  (3),  the 
power  is  found  to  be  32  kilowatts.   Details  of  the  theoretical 
power  available  calculations  are  shown  in  Appendix  A. 
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III.   EXPERIMENTAL  DEVELOPMENT 

A.   DESCRIPTION  OF  EXPERIMENT 

Reference  2  is  a  final  report  of  analytical  and  experi- 
mental programs  conducted  by  the  AVCO  Everett  Research 
Laboratory  to  determine  the  feasibility  of  using  a  cylindrical 
geometry  for  a  gas  dynamic  laser.   The  Naval  Weapons  Center, 
China  Lake,  California,  expressed  an  interest  in  the  concept 
of  a  cylindrical  GDL  and  requested  that  research  be  conducted 
to  determine  whether  a  cylindrical  shock  wave  could  be 
established  by  incorporating  the  recommendations  of  the  AVCO 
Everett  report.   Funding  of  this  project  was  provided  by 
Naval  Air  Systems  Command,  Code  320. 

Before  determining  the  actual  shape  of  the  nozzle  rings 
(see  figures  (8)  and  (9)),  it  was  first  necessary  to  calculate 
the  throat  height,  that  is,  the  spacing  between  each  nozzle 
ring  in  the  stack.   This  spacing  is  determined  by  the  mass 
flow  that  could  be  supplied  to  the  nozzle  stack,  which  is 
limited  by  the  feed  holes.   See  figure  (14) . 

Once  the  throat  height  was  determined,  the  exit  to  throat 
area  ratio,  and  thus  the  exit  Mach  number,  was  found  from  the 
geometry  of  the  nozzle  rings.   The  axisymmetric  method  of 
characteristics  found  in  Ref.  3  was  then  used  to  determine 
the  contour  of  the  nozzle  walls.   Once  the  equation  describing 
the  wall  shape  was  calculated,  the  design  was  sent  to  the 
machine  facility  at  NWC ,  China  Lake,  where  the  nozzle  rings 
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were  fabricated.   Figures  (8) ,  (9) ,  (10) ,  (12)  and  (14)  were 
provided  by  NWC ,  China  Lake. 

To  prevent  separation  of  the  boundary  layer  at  the  shock 
location,  suction  was  used  to  bleed  off  the  low  energy  flow 
in  the  lower  portion  of  the  boundary  layer.   Extensive 
analytical  investigation  of  supersonic  boundary  layer  bleed 
was  conducted  and  is  found  in  section  III-D  of  this  report. 

The  loss  in  total  pressure  across  a  normal  shock  of  the 
strength  anticipated  resulted  in  total  pressure  downstream 
of  the  shock  (using  a  reservoir  pressure  of  ten  atmospheres) 
to  be  less  than  one  atmosphere.   Ejector  nozzles  were 
necessary  to  lower  the  back  pressure  seen  by  the  flow  exiting 
the  nozzle  rings  and  thus  allow  them  to  start.   Design  of 
the  ejector  apparatus  is  found  in  section  III-E  of  this 
report. 

Figure  (7)  shows  the  assembled  experimental  apparatus 
with  dimensions.   The  design  allowed  for  testing  of  up  to 
eight  nozzles  at  the  design  Mach  number.   The  assembly  shown 
was  attached  to  the  free  jet  located  in  building  230  of  the 
Naval  Postgraduate  School.   Air  was  supplied  to  the  free  jet 
from  a  2400  cubic  foot  tank  capable  of  pressure  up  to  three 
hundred  pounds  per  square  inch.   The  moisture  was  removed 
from  the  air  before  it  entered  the  tank  by  means  of  a 
hydryer.   A  model  06231  Lammert  Industries  vacuum  pump  was 
used  to  evacuate  a  200  cubic  foot  tank  to  provide  suction 
for  boundary  layer  bleed. 
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6.5" 


FIGURE  (7)   CYLINDRICAL  NOZZLE  TEST  RIG 
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B.   NOZZLE  FEED  SYSTEM 

1.   Throat  Height 

It  was  found,  as  expected,  that  the  throat  height 

was  the  most  critical  dimension  of  the  entire  design.   A 

relatively  large  throat  height  is  desirable  in  that  the 

ratio  of  the  nozzle  exit  area  to  the  throat  area  is  kept 

small,  thus  keeping  the  exit  Mach  number  reasonably  low. 

This  would  have  the  advantages  of  ease  in  starting  and 

smaller  total  pressure  loss  across  the  shock. 

A  large  throat  height  also  requires  large  mass  flow 

requirements,  and  the  amount  of  air  that  can  be  fed  to  the 

nozzle  rings  is  limited  by  the  inlet  holes  in  the  center  of 

each  ring.   Figures  (10)  through  (14)  show  the  dimensions 

and  photographs  of  the  nozzle  rings.   The  area  of  the  feed 

holes  must  be  greater  than  the  combined  throat  areas  of  the 

nozzle  ring  stack  to  prevent  choking  of  the  flow  at  the  feed 

holes.   Thus  if  many  nozzle  rings  are  to  be  tested,  an 

extremely  small  throat  height  must  be  used.   This  would 

result  in  a  number  of  cor.pl ications : 

(i)   An  excessively  large  Aexit/Athroat /  causing 

high  supersonic  or  hypersonic  flow  at  the  nozzle 

exit  and  thereby  creating  a  problem  in  starting 

and  a  large  loss  in  total  pressure  across  the  shock. 

(ii)   Extreme  heat  concentrations  at  the  throat  (not  a 
problem  in  cold  f low  testing) . 

(iii)   Difficulty  in  fabricating  spacing  washers  of  uniform 
thickness,  causing  unacceptable  variations  in  the 
flow  due  to  changes  in  the  throat  height. 
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FIGURE  (13)   DIMENSIONS  OF  NOZZLE  RING  (END  VIEW) 
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The  throat  height  must  satisfy  the  continuity 
equation : 


m.   =  I     rn 
in   n=l  n 


where  n  is  the  number  of  nozzle  rings.   This  can  also  be 
expressed  as 


(2pinUinAin)+  =  6P*U*A*   •  d> 


Rewriting  equation  (1)  using  the  equation  of  state  and 
solving  for  the  throat  height: 


h  =  6¥MinAin(-p^)(f— > 

in 


Using  the  value  of  inlet  area  previously  obtained  and  using 
the  isentropic  relations  for  pressure  and  temperature  in 
equation  (2) ,  the  throat  height  in  inches  can  be  expressed 
as 


h  ■  -092  "-T  •  <3> 

(1  +  •2Mi„)J 


4- 

The  factor  of  two  in  equation  (1)  results  from  the 
nozzle  stack  being  fed  from  both  ends. 
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Equation  (3)  shows  that  the  throat  height  is  a  function  of 
inlet  Mach  number  only.   Details  of  the  throat  height 
calculations  are  found  in  Appendix  C. 

The  exit  Mach  number  is  determined  by  the  ratio  of 
the  exit  area  to  the  throat  area.   This  ratio  is  found  from 
the  geometry  of  figure  (11)  to  be 


Aexit    2^e(h+t) 


or 


2  7irth 


Ae    1.5(h+.5)  ... 

__        .  (4) 


Also,  assuming  isentropic  flow  through  the  nozzle,  the  area 
ratio  can  be  expressed  as 


A  Y+l 

_e  =  [  _L  [(  2  )( l+Izl  M  2n2TFTr  ,  (5) 

A    l  M   LVy+l      ^^  e  '  J        J  K    ' 


Equations  (4)  and  (5)  show  that  the  exit  Mach  number  is 
directly  dependent  upon  the  throat  height.   Since  the  throat 
height  is  a  function  only  of  the  inlet  Mach  number,  the  exit 
Mach  number  is  also  a  function  of  the  inlet  Mach  number. 
The  results  of  equations  (3) ,  (4)  and  (5)  for  various  inlet 
Mach  numbers  are  listed  in  Table  (1) .   Figures  (15)  through 
(17)  show  the  results  of  Table  (1)  in  graphical  form. 
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Table  1   Throat  Height  and  Exit  Mach  Number  for 
Various  Inlet  Mach  Numbers. 


Minlet  h  <inches>  Mexit 


.2  .0179  5.65 

.4  .0334  4.95 

.6  .0446  4.62 

.8  .0511  4.47 
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A  throat  height  of  .05  inches  was  chosen  to  allow 
for  the  lowest  exit  Mach  number  that  would  not  result  in  a 
choked  inlet  condition.   From  figure  (17)  ,  a  throat  height 
of  .05  inches  results  in  an  exit  Mach  number  of  4.5.   The 
axisymmetric  method  of  characteristics  was  used  to  design 
the  contour  of  the  nozzle  rings.   These  calculations  are 
found  in  section  III-C. 

Knowing  the  throat  height  between  each  nozzle  ring, 
the  mass  flow  rate  can  be  determined.   For  each  nozzle, 


m    -,      =   p*U*A* 
nozzle 


Using  the  equation  of  state  and  rearranging  yields 


P* 
Po(P~)  (49-1)     .2tt(1)  (105) 

A   =  [—2- -]   [ ]  . 

T*  h  144 

R(to  $-)  2 
o 

Substituing  the  appropriate  values,  the  mass  flow  per 
nozzle  is  found  to  be 


n^  =  1.08  lb  m/sec 
Total  mass  flow  is 

Atotal  =  6  (*n) 


m.  4.  i  =  6.4  8  lb  /sec 
total  m 
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2.   Axial  Variation  of  Mach  Number  and  Static  Pressure 
After  the  flow  enters  the  nozzle  stack,  changes  in 
Mach  number  and  static  pressure  occur  in  the  axial  direction. 
Shapiro  [Ref.  4]  develops  influence  coefficients  that  relate 
changes  in  one  flow  condition  to  changes  in  another. 
Assuming  constant  molecular  weight,  the  change  in  Mach 
number  is  related  to  the  area  change  by 


dM2      2(1+I^M2) 

I-7H     •  \v) 


vr  1-M2     A 


Appendix  D  shows  the  model  used  for  the  axial  flow  through 
the  nozzle  stack.   The  pressure  change  is  also  related  to 
changes  in  area  and  is  given  by 


dP  _  yM2   dA  fn, 

-p"  "  7— T  T     '  (7) 

1-M 


Using  the  geometry  of  the  nozzle  stack,  and  rearranging 
equations  (6)  and  (7) ,  axial  Mach  number  and  pressure 
changes  are  found  to  be 

1  +  l£  M  2 

M..,  =  M.  -  M.  t ,     (.06667)]       (8) 

l+l     1     l     i-M. 


and 


2 

P_.  =  P.  +  P.  [ i-y  (.06667)]   .         (9) 

1  +  1  1  1      1  \A       *> 

1  -  M . 
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A  computer  program  was  written  using  equations  (8)  and  (9) 
to  determine  values  of  Mach  number  and  pressure  in  the 
axial  direction.   Figures  (18)  and  (19)  are  the  results  of 
this  program  in  graphical  form.   A  listing  of  the  program 
and  the  results  are  found  in  Appendix  D,  along  with  details 
of  the  preceding  development. 

3 .   Radial  Variation  of  Pressure  and  Mach  Number 

As  the  flow  leaves  the  nozzles,  it  experiences 
increasing  area  as  it  moves  out  radially.   This  area 
increase  results  in  increasing  Mach  number  and  a  corresponding 
decrease  in  static  pressure. 

From  the  nozzle  ring  geometry, 


'       27Tr,  h 


A/A*  =  11.0  (r)   .  (10) 

Assuming  isentropic  flow: 

Y+l 

VA-i^d^M2)]^^  (11) 

and 

Y 
Po/P  =  [(l+I^M2)]^   .  (12) 
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Equations  (10)  ,  (11)  and  (12)  relate  Mach  number  and 
pressure  to  radial  position.   Table  (2)  is  a  summary  of 
these  three  equations  for  different  radii.   Figures  (20) 
and  (21)  are  graphical  representations  of  Mach  number  and 
pressure  as  functions  of  radius. 

C.   NOZZLE  RING  DESIGN 

1.   Nozzle  Contour 

As  was  shown  previously  in  Section  III-A,  an  exit 

Mach  number  of  4.5  was  chosen  for  design.   The  curvature  of 

the  nozzle  rings  was  determined  by  the  use  of  the  axisym- 

metric  method  of  characteristics. 

a.   Initial  Turning 

In  a  two-dimensional  nozzle  the  flow  is  turned 

through  (v,_.   ,/2)  degrees  at  the  throat,  where  v^ .   .is 
'    final      r  '        final 

the  value  of  the  Prandtl-Meyer  function  corresponding  to  the 
design  exit  Mach  number.   In  a  three-dimensional  case  the 
flow  expands  from  throat  to  tip  (and  therefore  the  Mach 
number  increases)  due  solely  to  the  increase  in  area.   Thus 
to  achieve  the  desired  exit  Mach  number  in  the  three  dimen- 
sional case,  the  flow  should  be  turned  less  than  in  the 
two-dimensional  case. 

Assuming  parallel  plates,  the  ratio  of  the  exit 
area  to  the  throat  area  reduces  to 


VAthroat  "  1-5'1 
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Table  2   Radial  Variation  of  Mach  Number  and  Pressure 


r  (in.)  M 


1.5  4.50 

2.0  4.84 

2.5  5.089 

3.0  5.239 

3.5  5.479 

4.0  5.674 

4.5  5.869 


Pstat    (P£ 

if) 

73, 

.11 

48, 

.30 

36, 

.07 

30, 

.39 

23, 

,27 

18, 

.85 

14, 

.92 
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Using  a  design  Mach  number  corresponding  to  an  area  ratio 
of 

A/A*  =  (A/A*)M  =  4.5(1/1.5) 
or 

A/A*  =10.4   , 

and  then  treating  the  initial  turning  as  if  it  were  a 
two-dimensional  case  should  give  approximate  results  (that 
is,  an  exit  Mach  number  near  4.5).   For 

A/A*  =10.4  , 


M   ..  =  3.96 
exit 


which  yields 


vexit  "  65-3° 


Therefore  the  initial  turning  of  the  flow  at  the  throat 
should  be 


v   . , 
A6  =  ex   =  32.65° 


where  6  is  measured  from  the  flow  direction  entering  the 
throat. 
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b.   Axisymmetric  Method  of  Characteristics 

Liepmann  and  Roshko  [Ref.  3]  develop  expressions 
for  finding  the  Prandtl-Meyer  function  (v)  and  the  flow  for 
direction  (0)  for  an  axially  symmetric  flow  and  are  given  by 


v3  =  ir(vl  +  v2>  +  l(ei-62' 


1        sin  9,  sin  6 

+  ^[sin  Ul  ~7—  A513  +  sin  v2   —?—   An^] 


(13) 


and 


63  "  J(vl-V2»  +  J(ei  +  e2' 


,        sin  81  sin  62 

+  ^[sin  yx  —F—   A?13  -  sin  v2   -~^—  An23  j 

(14) 


where  all  angles  and  subscripts  are  as  shown  in  figure  (22) . 

Using  equations  (13)  and  (14)  the  nonsimple  core 
of  the  characteristic  mesh  was  calculated.   Figure  (23)  shows 
the  entire  characteristic  mesh.   The  initial  turning  at  the 
sharp-edged  throat  was  divided  into  eight  segments  as  figure 
(23)  illustrates.   Section  II  of  Appendix  E  contains  tables 
of  Mach  number,  Mach  angle,  flow  angle  and  Prandtl-Meyer 
function  for  each  point  in  the  mesh. 

As  the  last  of  the  right  running  characteristics 
from  the  throat  reached  the  plane  of  symmetry  (point  70) , 
the  Mach  number  was  found  to  be  4.47,  and  the  flow  angle  was 
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x> 


+  0=90 


9  =  0' 

FIGURE  (22)   CHARACTERISTIC  ELEMENT 
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90°.   These  flow  conditions  show  that  by  reducing  the 
design  Mach  number  by  the  method  previously  described,  and 
then  treating  the  initial  turning  as  a  two-dimensional  case, 
accurate  results  will  be  obtained. 

c.   Establishing  Points  on  the  Nozzle  Wall 

The  left  running  characteristic  from  point  70 
is  a  straight  line,  and  the  flow  conditions  at  any  point  on 
this  characteristic  will  be  the  same  as  at  point  70,  that  is, 
M  =  4.47  and  0  =  90°.   Points  on  the  nozzle  wall  were  deter- 
mined by  calculating  points  along  a  right  running  character- 
istics in  the  -E,   direction  starting  from  a  point  on  the  left 
running  characteristic  between  points  70  and  76. 

Following  a  development  similar  to  Ref.  3, 
equations  (13)  and  (14)  were  modified  to  calculate  points 
of  the  characteristic  mesh  proceeding  in  the  -E,   direction. 
The  modified  equations  are 


v3  =  y(v1  +  v2)  +  j(e-L-e2) 


,         sin  61              sin  62 
+  T[-sin  y,  — A£,  ,  +  sin  y~  — — An„] 


1  — E, "^13     ll    M2    r0    ""23 

(15) 


and 


63  =  J(vl-V2}  +  y(6l  +  e2) 


,         sin  e  sin  62 

+  ^[-sin  Wl  — _  A^13  -  sin  y2  _—  Ar^l 


(16) 
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where  angles  and  subscripts  are  as  shown  in  figure  (E-l) 
found  in  Appendix  E.   Equations  (15)  and  (16)  were  used  to 
construct  right  running  characteristics  in  the  -£  direction 
from  points  72,  74,  76  and  78  (not  shown)  until  the  left 
running  characteristics  reflected  from  the  centerline  were 
met.   Starting  at  the  throat  the  wall  was  drawn  at  32.65° 
(relative  to  the  centerline)  until  point  19  was  reached, 
where  the  wall  angle  was  reduced  to  the  flow  angle  that  was 
calculated  for  point  19.   This  procedure  was  continued  for 
points  32,  37  and  50.   The  last  point  (point  80)  was  the 
result  of  changing  the  angle  of  the  wall  to  the  estimated 
flow  angle  at  that  point.   This  line  was  extended  until  a 
radius  of  1.5  inches  was  reached.   Limiting  the  exit  radius 
to  this  value  was  an  effort  to  keep  the  exit  Mach  number 
as  low  as  possible  as  was  discussed  previously, 
d.   Wall  Shape 

The  r  and  z  coordinates  of  the  points  where  the 
wall  slope  changes,  as  well  as  the  throat  and  the  exit  points, 
were  fit  with  a  fourth-degree  polynomial  using  the  Hewlett 
Packard  model  9830-A  computer  and  a  polynomial  regression 
program.   The  coordinates  of  the  points  and  the  coefficients 
of  the  fourth-degree  polynomial  are  shown  in  Table  (3) . 
Figure  (24)  shows  the  nozzle  contour  from  throat  to  exit. 
2.   Trailing  Edge 

As  can  be  seen  from  figure  (24)  ,  the  nozzle  contour 
is  not  tangent  to  the  ring  centerline  at  r  =  1.5  inches. 
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Table  3    Description  of  Wall  Contour 


A.   Wall  Point  Coordinates  from  Method  of  Characteristics 


Point  Number 

X  Coord. 

(throat) 

1.00 

19 

1.06 

32 

1.11 

37 

1.18 

50 

1.29 

80 

1.44 

(exit) 

1.50 

Y  Coord. 

.025 

.066 

.097 

.138 

.173 

.215 

.226 

B.   Coefficients  of  Fourth  Degree  Polynomial  Fit  to  Points 

B(0)  =  -6.375 

B(l)  =  18.144 

B(2)  =-19.773 

B(3)  =   9.939 

B(4)  =  -1.910 

R-Square  =   0.9999  7 
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Rather  than  change  the  shape  of  the  wall  to  a  contour  that 
would  result  in  a  sharp  trailing  edge,  it  was  decided  to 
maintain  the  design  contour.   This  resulted  in  the  flow 
exiting  the  nozzle  at  a  5°  angle  and  a  blunt  trailing  edge 
as  shown  in  figure  (25) . 
3.   Throat  Shape 

A  sharp-edged  throat  was  used  to  obtain  the  highest 
degree  of  turning  for  the  shortest  length.   Reference  5 
states  that  to  be  assured  of  uniform  flow  at  the  nozzle, 
entrance,  the  radius  of  curvature  of  the  inlet  should  be 
two  to  three  times  the  throat  height.   Figure  (26)  shows  the 
details  of  the  nozzle  inlet  and  throat.   Due  to  the  original 
geometry  of  the  nozzle  ring  (before  machining)  the  ratio  of 
the  radius  of  curvature  to  the  throat  height  was  limited  to 
between  1.5  and  2.0.   The  shaded  area  in  figure  (26) 
represents  the  machining  tolerance  of  the  inlet  contour. 
The  center  of  curvature  was  displaced  to  the  left  of  the 
sharp  corner  to  assure  horizontal  tangency  of  the  flow. 

D.   EJECTOR  NOZZLE  DESIGN 

For  a  reservoir  stagnation  pressure  of  ten  atmospheres 
and  Mach-five  flow,  the  total  pressure  downstream  of  the 
normal  shock  will  be  less  than  one  atmosphere.   To  lower 
the  back  pressure  and  thus  allow  the  cylindrical  nozzles  to 
start, ejector  nozzles  were  designed  and  installed  in  the 
exhaust  section  of  the  apparatus  as  illustrated  in  figure 
(7).   Using  the  continuity,  momentum  and  energy  equations, 
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FIGURE  (26)   DETAIL  OF  NOZZLE  RING  THROAT 
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solutions  to  the  flow  were  found  for  different  combinations 
of  exhaust  area  and  ejector  nozzle  diameter.   A  computer 
program  was  written  to  solve  the  three  equations,  and  the 
results  are  tabulated  in  Appendix  F.   Figure  (27)  is  a  graph 
of  these  results.   Any  combination  of  nozzle  diameter  (Dl) 
and  exhaust  width  (D3)  which  results  in  a  chamber  pressure 
(PI)  below  the  total  pressure  line  will  allow  the  cylindri- 
cal nozzles  to  start.   A  nozzle  diameter  of  1.75  inches  and 
exhaust  width  of  6.5  inches  were  chosen.   This  combination 
results  in  a  back  pressure  of  .562  atmospheres.   Complete 
details  of  this  development  are  found  in  Appendix  F. 

E.   BOUNDARY  LAYER  BLEED 

Boundary  layer  bleed  was  employed  to  keep  the  flow  from 
separating  at  the  shock  location.   By  bleeding  off  the  low 
energy  air  in  the  boundary  layer  the  remaining  flow  will 
have  sufficient  energy  to  support  the  pressure  rise  associ- 
ated with  a  shock  of  the  strength  anticipated.   A  search  of 
the  literature  relating  to  boundary  layer  failed  to  locate 
specific  information  on  supersonic  boundary  layer  bleed 
techniques,  specifically,  hole  size  and  pressure  ratio 
required  to  remove  a  certain  fraction  of  the  boundary  layer. 
Since  information  could  not  be  found,  a  model  was  developed 
to  make  such  a  prediction. 

1.   Flow  Model 

A  model  fluid-  discussed  by  Stewartson  [Ref .  7]  was 
used  in  this  development.   The  fluid  is  assumed  to  have  a 
Prandtl  number  of  unity  and  the  viscosity  is  proportional  to 
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the  static  temperature  so  that  yp  =  y.p.,  where  subscript  1 
denotes  the  reference  condition.   Assuming  the  fluid  has  a 
Blasius  profile  and  that  stagnation  temperature  is  constant 
through  the  boundary  layer,  the  incompressible  profile  was 
transformed  to  a  compressible  profile  by  the  method  described 
in  Ref.  7. 

For  analytical  ease,  a  Blasius  profile  was 
approximated  by 

u/U  =  sin(Ari) 
where  A  =  54°  and  n  is  given  by  Ref.  7  as 

n  =  xV2/Re 

The  reference  length  is  denoted  by  X.   The  model  shown  in 
figure  (28)  was  used  to  describe  the  flow.   By  applying  a 
bleed  plenum  pressure  (P2)  lower  than  the  freestream 
pressure  (P, ) ,  a  streamline  (4)  is  assumed  to  stagnate  on 
the  downstream  edge  of  a  slot  of  width  H.   All  flow  below 
this  streamline  goes  through  the  slot,  which  is  assumed 
long  enough  to  have  parallel  streamlines  at  the  slot  exit. 
By  assuming  isentropic  flow  in  each  streamtube  below  (4) , 
it  is  possible  to  calculate  flow  conditions  at  the  slot 
exit  by  specifying  (P,/P2) .   Using  the  continuity  equation 
leads  to  a  slot  width,  H.   The  complete  set  of  equations 
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programmed  for  the  Hewlett  Packard  model  98  30  computer, 
and  is  listed  in  Appendix  G. 

2 .  Results  of  Calculation 

Figure  (29)  illustrates  the  results  of  the  computer 
program  and  plots  the  fraction  of  the  boundary  layer  that 
enters  the  slot  (Y/6)  as  a  function  of  the  pressure  ratio 
(P,/P2) ,  with  nondimensional  slot  width  H/6  as  a  parameter. 
Table  (4)  gives  the  results  for  the  flow  conditions  antici- 
pated for  radial  flow  at  a  Mach  number  of  4.5.  The  symbols 
used  are  shown  in  figure  (2  8)  or  have  been  given  previously 

with  the  exception  of  a   ,  which  is  the  stagnation  speed 

e 
of  sound  at  the  edge  of  the  boundary  layer. 

The  dot-dash-dot  line  in  figure  (29)  represents  the 

locus  of  maximum  Y/6  for  a  given  slot  width  and  pressure 

ratio. 

3 .  Bleed  Hole  Design 

Using  the  values  from  table  (4) ,  (specifically  H  and 
X)  the  total  bleed  area  was  found  to  be 


A.  .   ,  =  2ir(r  +X)H 
bleed       e 


Dleed 


=  2tt(1.5  +  3)  (.00143) 
=  .38  in 


Three  hundred  holes  were  drilled  through  each  plate  at  a 
diameter  of  .141  inches  as  shown  in  figures  (30)  and  (31) 
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FIGURE  (29)   PLOT  OF  FRACTION  BOUNDARY  LAYER  BLEED,  Y/6, 
AS  A  FUNCTION  OF  PRESSURE  RATIO  AND  NON- 
DIMENSIONAL  BLEED  SLOT  WIDTH  (FJ  M  REF.  8) 
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Table  4   Flow  Conditions  for  Boundary  Layer  Bleed  (from  Ref .  8) 


M  =  4.5              Re  =  500,000  x  =  0.25  ft 

Y  =  1.4             A  =  54°  aQ  =  1157  ft/sec 

e 

Ue  =  2318  ft/sec     p1  =  200  psf  P2  =  50  psf 

6  =  0.00504  ft      Y/6  =  0.385  H/6  =  0.283 

H  =  .00143  ft       m.  .  ,  =  1.82  E-4  slugs/sec 

total 


STREAMTUBE 

NUMBER  n  Y  AY  AH  M  M^,  m  AH/AY 

ft3  ^  ^ Slug5/Sec 

0  0 

1  .40  .60  .076         1.56         0.74         0.15 
.16           .40 

2  .40  1.88  .227         1.58         2.21         0.47 
.32           .80 

3  .39  3.01  .381         1.63         3.67         0.77 
.48         1.19 

4  .38  3.98  .537         1.69         5.10         1.04 
.64         1.58 

5  .37  4.75  .698         1.78         6.51         1.28 
.80         1.95 
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FIGURE  (30)   BLEED  HOLE  DESIGN 
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The  total  area  of  the  holes  on  one  plate  is  larger  than 

the  calculated  bleed  area.   This  was  done  for  three  reasons: 

1)  to  bleed  off  a  greater  fraction  of  the  boundary 
layer  than  calculated  if  needed  (see  figure  (29)), 

2)  to  keep  the  same  bleed  area  but  vary  the  radial 
position  of  the  effective  bleed  area  (i.e.,  if 
bleed  was  desired  at  an  inner  radius,  the  outer 
holes  were  covered,  or  vice  versa.), 

3)  since  the  bleed  plenum  was  fitted  with  pressure 
taps,  flow  velocities  could  be  attained  by 
securing  the  suction  to  the  plenum  and  recording 
the  pressure  drop  in  the  plenum  due  to  the  flow 
over  the  bleed  holes.   By  covering  all  but  one 
row  of  holes  the  velocity  at  a  desired  radius 
could  be  obtained  by  using  the  isentropic  relation 
for  pressure  and  Mach  number  or  Bernoulli's 
equation. 

F.   EXPERIMENTAL  PROCEDURE 

Initially  six  nozzle  rings  were  tested  in  the  apparatus 
illustrated  in  figure  (7) .   Stagnation  pressure  applied  to 
the  nozzle  stack,  controlled  by  means  of  a  regulator  valve, 
was  increased  to  a  maximum  of  130  pounds  per  square  inch. 
Boundary  layer  bleed  was  attained  by  connecting  the  vacuum 
pump  directly  to  the  bleed  plenums.   The  flow  was  observed 
using  shadowgraph  and  Schlieren  techniques.   A  typical  run 
consisted  of  establishing  the  sensitivity  of  the  Schlieren 
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(or  shadowgraph)  system  to  the  desired  level,  starting  the 
vacuum  pump  and  allowing  it  to  reach  its  lowest  pressure, 
and  then  increasing  the  pressure  to  designated  levels  where 
photographs  were  taken. 

The  number  of  nozzles  was  then  changed  from  six  to  four 
and  the  spacing  of  the  nozzle  rings  (the  throat  height)  was 
increased  from  .05  inches  to  .09  inches.   This  change 
decreased  the  exit  Mach  number  from  4.5  to  3.9.   Since  the 
total  pressure  loss  across  the  shock  was  less,  the  back 
pressure  did  not  need  to  be  lowered  to  below  atmospheric. 
Therefore  the  ejector  nozzles  were  removed.   Removing  the 
ejectors  allowed  for  longer  runs  at  high  pressures  due  to 
a  decrease  of  approximately  forty  percent  in  the  mass  flow 
from  the  reservoir. 

The  four  nozzle  rings  were  tested  in  the  apparatus 
modified  as  shown  in  figure  (32)  in  a  fashion  similar  to  that 
previously  described  for  six  nozzles.   One  important  differ- 
ence was  the  use  of  a  200  cubic  foot  tank  to  increase  the 
bleed  capability  of  the  vacuum  pump.   The  pump  was  connected 
to  the  tank  and  was  able  to  lower  the  pressure  inside  it  to 
approximately  one-half  pound  per  square  inch.   The  bleed 
plenums  were  then  connected  to  the  bleed  tank  and  the  suction 
was  controlled  with  a  gate  valve. 

Flow  conditions  were  observed  and  photographed  at 
pressures  up  to  the  maximum  attainable  pressure  of  250  pounds 
per  square  inch.   Photographs  were  taken  with  and  without 
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boundary  layer  bleed  to  compare  the  resulting  flow  fields. 
The  bleed  hole  pattern  was  also  altered  to  obtain  pressures 
in  the  bleed  plenum  due  to  the  velocity  of  the  flow  across 
the  holes  at  various  radii  as  discussed  in  Section  III-E. 
A  pointed  probe  was  inserted  into  the  flow  at  various 
locations  and  different  total  pressures  to  determine  regions 
of  subsonic  and  supersonic  flow. 
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IV.   RESULTS  AND  CONCLUSIONS 

A.   RESULTS 

As  was  previously  mentioned,  six  nozzle  rings  were 
initially  installed  in  the  test  apparatus  shown  in  figure 
(32) .   Normal  shocks  were  obtained  in  all  except  the  end 
nozzles,  but  attempts  at  moving  the  shocks  beyond  the  nozzle 
exit  to  form  a  "barrel"  shock  were  unsuccessful  at  pressures 
of  up  to  130  psi.*   Figures  (33) (a)  and  (b)  are  Schlieren 
photographs  of  the  six  nozzles  at  pressures  of  70  and  130 
psi,  respectively.   Before  attempting  higher  pressures  it 
was  decided  to  increase  the  throat  area  in  order  to  lower 
the  exit  Mach  number  and  thus  make  it  easier  to  move  the 
shocks  out  of  the  nozzles.   Since  the  sum  of  the  throat  areas 
must  remain  smaller  than  the  feed  area  to  prevent  choking  at 
the  ends  of  the  stack,  the  number  of  nozzle  rings  was  reduced 
to  four  and  the  throat  height  increased  to  .09  inches.   This 
resulted  in  an  exit  Mach  number  of  3.9  and  increased  the  Mach 
number  at  the  feed  holes  to  .9. 

Three  distinct  flow  patterns  were  observed  in  the  flow 
over  the  range  of  pressures  up  to  a  maximum  of  250  psi. 
Figures  (34) ,  (35)  and  (38)  are  Schlieren  photographs 
showing  the  development  of  the  first  stage  flow  structure. 
Figure  (37)  is  a  schematic  of  the  flow  just  prior  to  its 


* 

Applied  pressures  are  gage  pressures  unless  otherwise 

specified. 


75 


(a)  70  PSI. 


(b)  130  PSI. 


FIGURE  (33) 


SCHLIEREN  PHOTOGRAPHS  OF  NORMAL  SHOCKS 
STANDING  IN  SIX  NOZZLES 
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(a)  NO  FLOYi 


(b)  50  PSI. 


FIGURE  (34) 


SCHLIEREN  PHOTOGRAPHS  OF  FOUR  NOZZLES 
BEFORE  "STARTING" 
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(a)  75  PSI 


(b)  100  PSI.  -  FIRST  SHOCK 
STRUCTURE  ESTABLISHED 


FIGURE  (35) 


SCHLIEREN  PHOTOGRAPHS  OF  CENTER  THREE 
NOZZLES  STARTING 
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(a)   Flow  at  Pressure  Less  Than  50  PSI. 


COANDA 
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FEED 
HOLE 
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/-**"  /    SLIP-LINE 


(b)   Resulting  Coanda  Jet  Attached  to  Nozzle  Wall 


FIGURE  (36)   SCHEMATIC  OF  COANDA  JET  FORMATION 
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FIGURE  (3  7)   SCHEMATIC  OF  FLOW  AT  A  TOTAL  PRESSURE 
OF  100  PSI. 
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FIGURE  (3  8) 


SCHLIEREN  PHOTOGRAPH  OF  FLOW  JUST  BEFORE 
TRANSITION  TO  SECOND  FLOW  STRUCTURE  - 
TOTAL  PRESSURE  125  PSI . 
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transition  to  the  second  stage  at  145  psi.   It  can  be 
readily  observed  that  failure  of  the  end  nozzles  (actually 
"half"  nozzles)  to  start  complicates  the  flow  structure  a 
great  deal.   Coanda  jets  appear  to  be  attached  to  end  nozzle 
walls  creating  a  "bubble"  of  relatively  slow  flow.   A  slip 
layer  results  between  this  slow  region  and  the  Coanda  jet 
and  appears  in  the  Schlieren  photographs  as  a  light  area. 
It  was  originally  considered  that  these  bright  regions  could 
be  expansion  waves  turning  the  flow  back  to  each  wall. 
Inserting  a  probe  into  this  region  as  shown  in  figure  (42-b) 
established  that  the  flow  was  in  fact  subsonic  since  it 
failed  to  produce  oblique  shock  waves  from  the  tip  of  the 
probe.   Pressure  readings  were  also  recorded  on  the  wall  at 
the  nozzle  exit  and  by  assuming  isentropic  flow  from  there 
to  ambient  pressures  (calculating  upstream) ,  the  Mach  number 
was  calculated  and  found  to  be  subsonic.   Table  (5)  lists 
the  pressure  data  taken  and  the  calculated  Mach  numbers. 
Figures  (39)  and  (40)  show  the  pressure  and  Mach  number  at 
the  wall  as  functions  of  total  pressure.   The  slip-line 
causes  a  bright  area  in  the  Schlieren  photographs  due  to  the 
higher  entropy  and  temperature  generated  from  friction, 
leading  to  an  area  of  lower  density. 

Careful  analysis  of  the  photographs  of  the  flow  up  to 
pressures  of  approximately  145  psi.  (e.g.  before  transition 
to  the  second  phase)  show  a  shock  wave  starting  normal  to 
the  wall  within  the  end  nozzles  and  then  curving  to  the 
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Table  5   Wall  Static  Pressure  at  Nozzle  Exit  and 
Calculated  Mach  Number. 


Total  Gage 
Pressure  (psi) 

Wall  Static 
Pressure  (psi) 

* 

Mach  Number 

at  Wall 

0 

12.94 

.43 

20 

12.42 

.50 

40 

11.59 

.59 

60 

11.0 

.66 

80 

10.51 

.71 

100 

9.72 

.79 

120 

9.28 

.84 

130 

- 

/ 

8.79 

-   / 

.89 

135 

- 

/ 

8.67 

-   / 

.90 

140 

8.56 

/ 

10.63 

.91  / 

.70 

145 

10.66 

/ 

10.73 

.69  / 

.69 

150 

10.69 

/ 

10.78 

.69  / 

.68 

155 

10.79 

.68 

160 

10.78 

.68 

* 

The  second  values  listed  resulted  from  decreasing  the  total 

pressure  from  150  psi  (illustrating  the  hysteresis  effect) . 
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angle  of  the  nozzle  wall,  where  the  slip  layer  forms 
between  the  Coanda  jet  and  the  subsonic  "bubble."   The  slip 
layer  continues  to  curve  back  toward  the  end  wall.   Observable 
in  the  lower  jet  are  two  small,  nearly  normal  shock  waves  — 
one  where  the  shock  attached  to  the  wall  begins  to  curve 
outward  and  one  from  the  trailing  edge  of  the  nozzle  ring. 
These  combine  to  slow  the  flow  to  subsonic  velocities  as 
testing  with  the  probe  revealed.   The  jet  on  the  upper 
nozzle  is  smaller  and  faster  which  is  attributable  to  a 
slightly  smaller  opening  at  the  throat.   Figure  (38)  shows 
the  alternately  light  and  dark  regions  of  a  series  of  shock 
diamonds  in  the  upper  Coanda  jet,  indicating  a  higher  Mach 
number  than  in  the  lower  nozzle.   The  slip  line  in  the  upper 
case  also  bends  back  toward  the  wall. 

The  high  velocity  jets  occur  for  one  of  two  reasons. 
First,  since  the  end  nozzles  are  really  only  "half-nozzles," 
the  throats  are  only  half  the  thickness  of  those  in  the 
center,  in  this  case  on  the  order  of  .045  inches.   With 
these  small  openings  located  at  the  ends  of  the  nozzle  stack 
where  the  flow  is  entering  axially  at  near  sonic  conditions, 
it  has  a  difficult  time  "turning  the  corner"  to  exit  out  the 
first  throat  opening.   If  this  were  the  case,  whatever  flow 
did  make  the  90°  turn  may  attach  itself  to  the  nozzle  wall 
as  shown  below. 


vn-rn     \nzzz2ZZZ2 


86 


A  second,  more  feasible  explanation  is  that  regardless  of 
whether  or  not  the  flow  has  trouble  turning  the  corner,  the 
flow  that  travels  through  the  throat  starts  as  a  jet  of  air 
near  the  wall  at  low  pressures  (in  the  order  of  50  psi.)  is 
shown  in  figure  (36-a) .   This  jet  entrains  the  air  near  the 
nozzle  wall  and  thus  lowers  the  pressure  in  the  wall  region. 
The  jet  of  air  then  moves  toward  the  wall  and  attaches  to 
it,  displaying  the  "Coanda  effect."   As  the  flow  in  the  jet 
moves  out  radially  it  slows  and  increases  in  pressure, 
forcing  the  slip  line  back  toward  the  wall  and  creating  a 
"bubble,"  observable  in  figure  (35).   As  the  total  pressure 
increases  toward  the  transition  point  of  145  psi. ,  the 
bubble  becomes  longer  as  the  flow  in  the  jet  takes  a  longer 
distance  to  reach  a  higher  pressure. 

The  three  nozzles  in  the  center  established  normal 
shocks  quite  readily  and  moved  outward  to  the  nozzle  exit 
with  increasing  pressure.   As  can  be  seen  in  the  Schlieren 
photograph  of  figure  (3  8) ,  connecting  these  three  waves 
would  result  in  a  cylindrical  shock  bowed  inward  at  the 
center.   Reference  2  states  that  a  loss  of  total  pressure 
in  an  axial  direction  will  manifest  itself  in  a  shock  wave 
that  is  conical,  tapering  inward  in  the  direction  of 
decreasing  total  pressure,  rather  than  strictly  cylindrical. 
For  a  stack  of  nozzle  rings  fed  from  both  ends,  as  was  the 
case  in  this  experiment,  such  a  loss  in  total  pressure  would 
result  in  the  condition  observed  at  pressures  between  90  and 
140  psi.,  that  is,  a  tapering  inward,  from  both  ends.   Total 
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pressures  different  from  the  region  mentioned,  however,  did 
not  result  in  a  distortion  of  this  type  (see  figures  (41) 
through  (48)).   Seemingly,  if  axial  total  pressure  losses 
were  observable  in  a  concave  shock  at  pressures  near  100  psi. 
they  would  also  be  observable  at  much  higher  pressures 
rather  than  disappear  near  145  psi. 

At  a  pressure  of  approximately  ten  atmospheres,  transition 
to  the  second  phase  flow  structure  occurred.   Figures  (41) (a) 
and  (b)  are  Schlieren  photographs  of  the  flow  at  total 
pressures  of  150  and  175  psi.   This  transition  consisted  of 
a  sudden  rise  in  the  wall  pressure,  and  therefore  a  drop  in 
the  Mach  number,  as  shown  in  figures  (39)  and  (40)  ,  respec- 
tively.  The  Coanda  jets  no  longer  form  a  subsonic  bubble, 
but  instead  the  slip  layer  extends  upward  at  an  angle  of 
19  degrees.   The  cause  of  the  transition  to  the  condition 
shown  in  figure  (43)  is  the  movement  of  the  normal  shock 
segments  beyond  the  trailing  edges  of  the  nozzles.   Since 
the  flow  leaves  the  nozzles  at  approximately  a  five  to  ten 
degree  angle,  oblique  waves  bend  the  ends  of  the  normal 
shock  segments  outward  giving  them  a  concave  shape.   Shock 
diamonds  form  behind  the  oblique  waves  serving  to  return 
the  flow  to  subsonic  conditions.   AVCO  Everett  in  Ref.  2 
predicts  that  disturbances  traveling  downstream  from  the 
trailing  edges  would  only  affect  the  cylindrical  shock  wave 
locally,  causing  outward  protrusions  on  the  surface  of  the 
shock.   Though  the  oblique  shocks  are  rather  large 
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(a)  150  PSI 


(b)  175  PSI. 


FIGURE  (41) 


SCHLIEREN  PHOTOGRAPHS  SHOWING  CURVED  SHOCK 
SEGMENTS  BEYOND  TRAILING  EDGES  OF  NOZZLES 
AFTER  TRANSITION  TO  SECOND  PHASE 


39 


(a)  150  PSI.  -  SUBSONIC 
BEHIND  CURVED  SHOCK 
SEGMENT 


(b)  150  PSI.  -  SUBSONIC 
IN  WALL  REGION 


FIGURE  (42) 


SCHLIEREN  PHOTOGRAPHS  WITH  PROBE  SHOWING 
REGIONS  OF  SUBSONIC  FLOW 
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FIGURE  (43)   SCHEMATIC  OF  THE  FLOW  STRUCTURE  AT  A 
TOTAL  PRESSURE  OF  175  PSI. 
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FIGURE  (44) 


SCHLIEREN  PHOTOGRAPHS  OF  FLOW  JUST  BEFORE 
DEVELOPING  INTO  THIRD  STAGE  -  TOTAL 
PRESSURE  OF  200  PSI. 
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(a)  220  PSI 


(b)  240  PSI. 


FIGURE  (45) 


SCHLIEREN  PHOTOGRAPHS  SHOWING  DEVELOPMENT 
OF  FLOW  TO  THIRD  STAGE 
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"disturbances,"  this  prediction  appears  to  hold  by  the  fact 
that  the  ends  of  the  normal  shock  segments  "protrude" 
outward.   Flow  exiting  in  a  more  nearly  horizontal  direction 
would  minimize  the  strength  of  these  oblique  waves,  thereby 
making  the  shock  segments  less  curved  and  closer  to  joining 
into  a  complete  cylindrical  shock  wave. 

The  second  phase  is  also  characterized  by  the  disappear- 
ance of  the  subsonic  "bubble"  in  each  end  nozzle.   As  the 
normal  shock  in  the  passage  above  (or  below)  the  end  nozzle 
moves  beyond  the  trailing  edge,  a  pressure  mismatch  occurs 
between  the  very  low  pressure  region  behind  the  shock  and 
the  Coanda  jet  coming  off  the  trailing  edge  on  the  opposite 
side  of  the  nozzle  ring.   This  causes  the  Coanda  jet  to  be 
diverted  upward  (or  downward) ,  and  the  resulting  slip-line 
to  extend  upward  (or  downward) .   The  angle  that  the  slip-line 
forms  increases  from  19  degrees  at  transition  (150  psi.)  to 
29  degrees  at  the  maximum  pressure  of  250  psi.   This  change 
is  due  to  the  increase  in  the  pressure  mismatch  created  at 
the  trailing  edge.   Since  the  trailing  edge  is  blunt,  the 
Coanda  effect  will  continue  to  force  the  jet  to  higher  angles 
as  it  tries  to  move  around  the  trailing  edge.   This  is  not 
the  case  with  a  sharp  trailing  edge  where  the  angle  of  the 
Coanda  jet  would  remain  the  same. 

As  the  pressure  was  increased  from  150  psi.  to  200  psi. 
the  flow  maintained  the  basic  configuration  shown  in  figure 
(43).   Figures  (41),  (42)  and  (44)  are  Schlieren  photographs 


94 


of  the  flow  in  this  range  of  pressures.   At  200  psi  the 
oblique  waves  from  the  trailing  edges  are  longer,  better 
defined,  and  at  a  smaller  angle.   The  curved  shock  segments 
have  moved  to  approximately  a  quarter  of  an  inch  beyond  the 
trailing  edges.   At  approximately  210  psi.  the  flow  structure 
is  altered  a  second  time. 

Transition  to  the  third  stage,  unlike  the  abrupt  change 
at  145  psi.  ,  was  a  gradual  development  of  normal  shock  segments 
in  the  nozzle  wakes.   Between  pressures  of  205  and  215  psi. 
the  normal  segments  slowly  appear,  along  with  two  weak  oblique 
shocks  at  the  trailing  edges  where  one  strong  wave  had  been 
at  lower  pressures.   Increasing  the  total  pressure  to  250  psi. 
changed  the  flow  structure  only  by  moving  both  the  normal  and 
curved  segments  of  the  shock  slightly  farther  out  from  the 
trailing  edges.   Figure  (47)  is  a  schematic  of  the  fully 
developed  flow.   Figure  (46)  is  a  Schlieren  photograph  of 
the  flow  at  the  maximum  pressure  applied,  and  shows  that  the 
normal  portions  of  the  shock  are  at  a  distance  of  approximately 
three  eighths  to  a  half  of  an  inch.   This  distance  is  short 
of  creating  a  lasing  cavity  wide  enough  to  allow  the  theoret- 
ically available  power  in  the  flow  to  be  extracted,  which  for 
four  nozzles  is  1.7  inches.   For  a  spacing  of  the  nozzles 
s  ,   (s)/(r  -r  )  =0.85  ,   where  r   is  the  shock  location. 
As  a  fraction  of  the  minimum  shock  location  allowable  to 
extract  the  power  available,  (rm~re)  t    the  location  of  the 
shock  wave  is: 
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FIGURE  (46) 


SCHLIEREN  PHOTOGRAPH  OF  FLOW  AT  MAXIMUM 
TOTAL  PPESSURE  OF  250  PSI.   NORMAL  SHOCK 
SEGMENTS  APPROXIMATELY  3/8  INCH  FROM 
TRAILING  EDGES 
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FIGURE  (47)   SCHEMATIC  OF  THE  FLOW  STRUCTURE  AT  A 
TOTAL  PRESSURE  OF  250  PSI. 
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Boundary  layer  bleed  was  applied  over  the  range  of 
pressures  with  little  if  any  effect  on  the  flow.   Figures 
(48) (a)  and  (b)  are  Schlieren  photographs  of  the  flow  at 
200  and  250  psi.  respectively  with  a  one-half  psi. bleed 
pressure  applied  to  the  first  three  rows  of  bleed  holes  on 
the  upper  plate.   At  a  total  pressure  of  250  psi  the  flow 
entrained  between  the  slip-lines  is  diverted  toward  the 
bleed  plate.   This  was  the  only  noticeable  change  caused 
by  the  bleed  over  the  entire  range  of  pressures.   Attempts 
were  made  at  applying  enough  suction  to  start  the  end 
nozzles  but  were  unsuccessful. 

Homogeneity  of  the  flow  in  the  lasing  cavity  is  a 
necessity  for  an  operational  laser.   This  was  not  the  case 
for  the  fully  developed  flow  at  250  psi.  since  the  end 
nozzles  failed  to  accelerate  the  flow  to  supersonic  condi- 
tions.  Other  disturbances  such  as  the  oblique  shock  waves 
from  the  trailing  edges  should  not  be  present  to  produce 
homogeneous  conditions  in  the  lasing  cavity. 
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(b)  250  PSI. 


FIGURE  (48) 


SCHLIEREN  PHOTOGRAPHS  SHOWING  FLOW  WITH 
BOUNDARY  LAYER  BLEED  APPLIED  TO  UPPER  PLATE 
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B.   CONCLUSIONS  AND  RECOMMENDATIONS 

Testing  of  the  cylindrical  nozzle  stack  at  a  total 
pressure  of  250  psi,  established  a  shock  wave  approximately 
a  half  of  an  inch  from  the  nozzle  exit.   The  shock  wave 
consisted  of  two  segments  normal  to  the  flow  (actually 
shock  "bands"  around  the  nozzle  stack)  connected  by  shock 
segments  curved  inward.   Failure  of  the  end  "half-nozzles" 
to  start  resulted  in  subsonic  flow  near  the  end  walls  for 
the  entire  range  of  total  pressures  investigated.   This 
prevented  the  normal  shock  from  reaching  the  end  walls. 

Since  subsonic  flow  existed  at  the  end  walls  and  the 
shock  did  not  reach  the  walls,  the  use  of  boundary  layer 
bleed  had  very  little  visible  effect  on  the  flow  structure. 

Flow  visualization  became  very  distorted  at  high 
pressures  due  to  the  extreme  density  gradients  and  the 
cylindrical  geometry  of  the  shock  waves. 

The  following  recommendations  are  made  for  further 
research  on  this  subject: 

1.  Increase  the  diameter  of  the  nozzle  rings  to  allow 
for  more  mass  flow  to  enter  the  nozzle  stack  without 
decreasing  structural  integrity  and  without  choking 
in  the  feed  passages. 

2.  Decrease  the  thickness  of  the  nozzle  rings  to  lower 
the  expansion  ratio  and  thus  lower  the  exit  Mach 
number. 

3.  Ensure  sharp  trailing  edges  on  the  nozzle  rings. 
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4.  To  help  avoid  starting  problems  in  the  end  nozzles, 
use  half  of  a  nozzle  ring  (sliced  along  the  edge) 

at  each  end.   This  would  eliminate  the  "half-nozzle" 
arrangement  and  the  Coanda  jets. 

5.  Investigate  total  pressures  in  excess  of  250  psi, 

6.  Install  pressure  taps  on  the  center  nozzle  wall  to 
determine  Mach  number.   Optical  distortion  due  to 
the  large  density  gradients  and  the  cylindrical 
geometry  results  in  difficulty  observing  shock 
angles  from  a  probe  inserted  in  the  flow. 
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APPENDIX  A 
THEORETICAL  LASING  POWER  AVAILABLE 

Intensity  is  defined  as  radiation  power  per  unit  area, 
with  units  of  kilowatts  per  square  foot.   Assuming  a 
Gaussian  intensity  distribution  on  each  mirror  as  shown 
in  figure  (A-l)  the  intensity  is 


2 
I(r)  =  I  e"r  c 


To  satisfy  the  boundary  condition  of 


I  =  Ip  (1/e2) 


at  r  =  r   ,   then 


o  =  2/re2 


Therefore: 


-2r2/r  2 
Kr)  =  I  e  .  (A-l) 


Equation  (A-l)  yields  a  Gaussian  peak  centered  at  r=0.   To 
obtain  a  ring  with  Gaussian  cross  section  having  the  peak 
intensity  at  r=r,  define 

r  =  r  -  r 
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FIGURE  (A-l)   LASING  MIRROR  INTENSITY  LOADING 
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Inserting  this  into  equation  (A-l)  yields 


I(r)  =  I  e   re~r       kwatts/cm 


To  obtain  the  total  power  output,  the  intensity  must  be 
integrated  over  the  mirror  area. 


rm   2 
m  —   0 

Power  =      f      f  I   exp  [2-JL-Y2  ]   r  dr  d 

r   r.  r  -r 

e   0  e 


m       r-r  r—    2 
P  =  2  I    /   r  expt^-i-V2  ]   dr  .       (A-2) 

e  e 


An  alternate  method  of  finding  the  power  available  from  the 
flow  is 


^       /     jtt     i.    \     ,specific> 
Power  =  (mass  flow  rate)  (  c  ) 

energy 


Power  =  E  (kjoules/lb  )  m  (lb  /sec) 


P  =  E  m  kwatts 


or 


P  =  2«aLp  a0(pe/p0)(ae/a0)MeB   .  (A-3) 
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Equating  equations    (A-2)    and    (A-3) : 


rm 


Iw     /        r  expt^-VT]2   dr  =  roLpoa^  (po/pJ  (ao/an)MeE    , 

(A-4) 


p         ^    —  ¥  e  o  o  "e'  ro '   e'  o  e 

r«         r«~r 
e  e 


where  the  only  unknown  in  equation  (A-4)  is  the  outer  radius 
of  the  mirror,  r  . 

In  solving  the  integral  Q  on  the  left  side  of  equation 
(A-4)  define 


r-r 
x  =  — — 


re-r 


■V2~  , 


making 


dx  =  — dr 


and 


re-r 


r  -r 
r  =  e_  x  +  r 


Y2" 

Making  the  corresponding  changes  in  the  limits  of 
integration: 

@   r  =  re  ,   x  =  VT 

@   r  =  rm  ,   x  =  -V^   , 
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the  integral  now  takes  the  form 


-V?       r-r       2   r  -r 
Q  -    /    (r  +  -e—  x)  e~x   (-4=-)  dx  , 

VT  V^  V^ 


or 


-V2~  -     o          -    -V2~ 

v  r  -r     2       r  -r     v 

/  r  (-4=-)  e~x  dx  +  (-£=. -  )         f 

yr  V^                V^    yr            (A_5) 


/      (-4=-)  e~X  dx  +  (-1= -  )    /   x  e~x  dx  * 


Since  the  second  term  on  the  right  side  of  equation  (A-5) 
contains  an  integral  of  an  odd  function  between  symmetric 
limits,  the  integral  goes  to  zero  leaving 

V2"     (r  -r)      2 
Q  =  -2   /   F    e     e"x  dx   .  (A-6) 

0  yr 

Multiplying  and  dividing  the  right  side  of  equation  (A-6)  by 
it  /2   and  rearranging, 

Y2~7r  (r-r  )    V2    9      2 
q  = ®_    /   _£_  e  x   dx   .  (A-7) 

0   yv 

Noting  that  the  integral  in  equation  (A-7)  is  the  definition 
of  the  error  function  gives 

Q  ^yrr/2   r  (r-r  )  erf(/T)   .  (A-8) 
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Substituting   equation    (A-8)    for   the   integral    in  equation 
(A-4)  : 


Iph/viTr    (r-re)    erf(/T)]    =   reLpoao  (Pe/PQ)  (ae/aQ)MeE    , 


or,    solving   for  r    , 


r-    (re/2)    +    1/2  -^re2  +  C  (A-9) 


where 


4    r      L  ME 

c  =  — 2 2_- 

^"72   I     erf{\J2) 


Since 


r    +  r 
—         em 
r  =  * — 


equation    (A-9)    can  be  written 


r      =        r   2  +C         .  (A-10) 

m  e 


The  values  in  the  expression  for  C  are 


r   =1.5  inches 
e 

Po  =  (.002378)  (32.174)  (10)  =  .765  lbm/ft3 
aQ  =  1119  ft/sec 
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<a/ao>M=4.5    =    «445 
(^o>M=4.5    =    -01745 


NI      =    4.5 
e 


E  =  5  kwatts/(lb  /sec) 


1=3  kwatts/cm2  =  2.8  xlO3  kwatts/ft2 
erf(/T)  =  2(. 9773  -  .500)  =  .9546   . 


Using  these  values  in  the  expression  for  C: 


C  =  .2692 (L)  inches2 


For   r   =1.5  inches,  equation  (A-10)  becomes 


r   -  V2.25  +  .2692  L 


m 


where  L  is  the  length  of  the  nozzle  stack  in  inches. 
Dimension  L  is  illustrated  in  figure  (7) . 

As  an  example  the  power  available  from  six  nozzles  was 
calculated.   The  theoretical  lasing  power  available  was 

found  to   be 

r 

m  r-r       /—     2 

Power   =   2ttI  /  r  exp  [        _y2    ]'   dr 

P      r  r-r 

e  e 

=   2ttI      erf  (\[2    )    y]v/2       [r  (r-r    )  ] 
P  e 

=    2tt(3)  (30.52/144)  (.9546)  (tt/2)  **  [1 .  635  ( .  135)  ] 
Power   =    32.16    kilowatts. 
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APPENDIX  B 
CALCULATIONS  FOR  DETERMINING  INLET  (FEED)  AREA 

Figure  (13)  shows  the  dimensions  of  the  nozzle  ring 
The  area  of  the  hole  segment  between  two  centers  is 

A  =  length  x  width 

circumference    ,,  ,   ,.    .   . 
=  -« x  (hole  diameter) 

=  27T(3/8  +  3/4)  (1/8)  x  (3/8) 

A  =  .1407  in2   . 

Since  there  are  four  segments: 


A1   =  4(.1407) 

2 

A,  =  .662  7  in 


The  area  of  the  remaining  half  circles  is 


A2  =  a  x  jno^i, 

A2  =  .4418  in2  . 


Therefore  the  total  inlet  area  is 


Ain  :  Al  +  A2 


2 

A.   -  1 . 104  in 
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APPENDIX  C 
DETAILS  OF  THROAT  HEIGHT  CALCULATIONS 

The  feed  flow  and  the  flow  exiting  through  the  nozzles 
must  satifsy  the  continuity  equation: 


min  "   Zn  Nozzle 
n=i 


2PinUinAin  =  6P*U*A*  (C"1) 


where  the  factor  2  in  equation  (C-l)  results  from  the 
nozzle  stack  being  fed  from  both  ends. 

Using  the  equation  of  state,  and  noting  that  velocity 
is  the  product  of  Mach  number  and  the  speed  of  sound, 
equation  (C-l)  can  be  written 


P 

(M,  a, J(A.  )  =  3(J£.)(l)(a*)(2Trr.h) 


RT.     in  in'   in'     xrtx'  '  t 

in 


Pin  P* 

M.  A.   =  6  — h 

in  in 


/— —   in  in      /— r 
V*in  VT* 

Solving  for  the  throat  height 


in 
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Po  To  k 

Multiplying  the  right  side  by   (-— )  and   (;=— )  2    , 

o         o 

equation    (C-2)    can  be  written   as 


h=    (irir>    "in'^'^'f1^    (^S-)!'2        •  <c"3> 

o  o  in 


Using  the  isentropic  equation  relating  stagnation  temperature 
and  pressure  to  static  temperature  and  presure  respectively, 
and  applying  them  to  the  Mach  equals  one  condition, 
equation  (C-3)  becomes 


h=  <W»  Min(^)(37fT)(-8333)('r2-),S 
o  in 


or 


M.  (P.  /P  ) 
h  =  .0917  — — °    .  (C-4) 


^ii/V 


"2 


Noting  that 


!°  =  (!?)YA-1.  (1  +  i2-ilM2,Y/Y-l 


equation  (C-4)  reduces  to 


h  =  .0917  M.   -^-i —-t- 

1   (T  /T.  )   2  (T  /T.  ) T/T  X 
v  o'  in'    v  o'  in' 
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or 


Y+l 


=  .0917  M.  (T.  /T  )i(Y  1) 


in   in'  o 


-(Y+l) 

h  =  .0917  M.  (1  +  (Y:1)  M.  2,217=17 
in       2    in 


.0917  M 


h  = 


in 


(1+  .2  Min2)3 
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APPENDIX  D 
AXIAL  VARIATION  OF  MACH  NUM3ER  AND  PRESSURE 

From  Ref.  4,  influence  coefficients  for  flow  in  a 
channel  are  given  as 


dM2        (1  +  ^  r2)    dA 

^V  =  -  2  ^ °*  (D-l) 


M  1-M 


and 


dP    yM2  dA 

P  "  7~T  A 
1-M 


(D-2) 


assuming  isentropic  flow  and  constant  specific  heat  and 
molecular  weight.   Using  the  model  shown  in  figure  (D-l) 
and  an  inlet  area  of  1.105  square  inches, 


Aout    1 

"s:      "  IT 
in 


Aout  -  <*T>  (Ain) 


A   .  =  .06906  in2 
out 


Since 


A2   Al   Aout   ' 
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A2  "  Al  _     .0  6906 

~iq   *'  1.165  -  .OG906" 


or 


$£■   =  .06667 
A 


Substituting  this  value  into  equations  (D-l)  and  (D-2)  and 
rearranging  yields 


+  I^M.2 


1  +  ^   Mi 
M.^n  =  M.  -  M.  [ 0     (.06667)] 


1  -  M. 

1 


and 

2 


(D-3) 


YM. 

pi+i  ■  pi  +  pi  Ittt  (-06667)]'  <D-4> 

1  -  M. 

1 


A  program  was  written  for  the  Hewlett  Packard  model  9830-A 
computer  to  calculate  the  axial  variation  of  Mach  number  and 
pressure  using  equations  (D-3)  and  (D-4) .   The  program 
determines  Mach  number  and  pressure  at  stations  one  through 
thirteen  using  values  of  inlet  Mach  number  and  total 
pressure  (PSF)  as  inputs.   A  listing  of  the  program  AXIAL 
and  the  results  of  the  program  are  found  in  this  appendix. 


Table  D-l   Symbols  for  Program  AXIAL 


Computer  Notation       Algebraic         Description 


Ml  M.  Inlet  Mach  number 

in 

P0  P  Stagnation 

pressure  (PSF) 

M(  )  Mach  number  at 

subscripted 
station 

P(  )  Static  pressure 

at  subscripted 
station 


116 


LISTING  OF  COMPUTER  PROGRAM  "AXIAL" 


i-n    Mi 


"::'  0     !: 


1.  0 

i  i  e 

1 2 


.1    -!• 

■I  !■:■; 
1  .:::: 


RI  NT    "  I  HI.  E 
R I  NT- 
INT 
PR  INI 
::,RIN 

:  j-  i  p 
j  ■••■  -i  ... 


p  i'  : 

i  E!  '/   \ 

•(■.IT! 


::  J 


i    i 


b  i  !•••!    ri-     n     I'liiL-rvi-i 


:"!  .• 


vl  I" 


j  ?  r  i  i  J 


■3  b  b  t. ."  * 


.i. 


■  i  i  ..    |.-'l  I-,  !- 


,-    .1 


'fclj      h 
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RESULTS  OF  PROGRAM  "AXIAL" 


:    ILET    MRCIv 


'  ■  1   ■!  C 


1  f    .  1 


MflCH# 


i 


0  „  1 

■  ■  .  :.    : 


■  '■) 


i  y 

3  2 


:  ..  .  ,  '  '..:' : '  ' 9  •.::'  £ 
U  ,.  I  3  o  1  o  A  I  . :! 
■'  i  :.;:;i  46287 
'..  .  6  98321 

-  „  <-    i  j  c  o  :'  '■■■  ■':.  '■■ 


21851] 6353  ; 
21868.  3449  3. 


!  !  1  I....  !....  !     Mi'1..!!': 


\ 


i  !'!    It 


STflTIC<PSF) 


4 
5 


i  •;■■■ 
1  •-■ 


8  ,  ~'f 

0      67286'   3  ! 

8,  23880  j  ■  ' 

8^ 1985  8  i 
8= 134687414 
0. 171 852 51 A 

0  .  j  '  ;;i  9  9  7  t ,  7  0  0 


[8951 .     : 
»56S. 77836 

i  988*1  =  6;;  4i-::;  i 


££3S'4. 

16139 

2  y  \  ■ 

341  48 

283 1 8 . 

;::,  ,  !  j 

i  '  4  2 . 

■  2  !.  1 

£0549. 

339 1  ': 

i  19;  . 

28728 » 

i  ';:  i 

c   i!  47=  55b1 18 
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RESULTS  OF  PROGRAM  "AXIAL"  (CON'T) 


!  i  i  ; 


•i  i 


H     # 


MRl>!# 


P    STfi  i'lC  !'  PSF 


1J   »    b 

,-.   ^   ,;::;  3  g  Q  Q.  |     g  5 

■I   r'"i  i"j  ~'   .H   i~-l  ■"'   ••'  |::" 

9  ,.  -1  3695' ' 


186     '        '462 
0     ! 


;  •  h 


bi 


3  y 

i  1 


9 ,,  £  4  7  5  9  6 1 6  4 


:jf  lb 

■  3  '  i : 

£  0  5 1  8  »  6  5  9  0  6 


N|  E!    MflCHtt 


!  U  I  ML     !■  Kt  :::•  b~"     d  '  I  bfc 


'IHl  ,HS 


P    8  T  fl  T  I  C  ( P  8  F 


1  ■'.-■ 
1 3 


;- '  -.  b-  ■  2  'o  tf  y  D  -:i  -J 
8. 556855 

9.4536 

P  ,  4 1 3  5  I :  ■  ■ 
s7.  379!  ■:•  1 "" 
8, 348754325 
0. 321 639263 
0  ..  2  9  7  £  2  6  3  3  9 
?, 275185     ■ 

Ci      -Ti  "J.  £•  cr  ..i  cr  q  ,•:  q 


17194, 1841 
17915.50  ' 
18472.03545 
1891 7  7407 
1 9  £  8 1  =  6  i  i; 

1  98  • :;  6 ,  8  2  3 ::;:'  2 
26S58, 478 1 
28271 ,82379 

2  0  3  8  6  =  4  A  • 

28518,74133 
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APPENDIX  E 
DETERMINING  POINTS  ON  THE  NOZZLE  WALL 

Following  a  development  similar  to  Ref.  4,  but  changing 
the  characteristic  network  as  shown  in  figure  (E-l) 
equations  can  be  developed  to  determine  points  along  a 
right  running  characteristic  proceeding  in  the  -£  direction. 

Along  the  characteristic  segment 

3  3 

/  |_  (v-e)  dn  -  /  sin  y  *2£L±  dn  (E-i) 


2 


and 


3  a  3 

/  ^  (v-9)  d£  =  -  /  sin  u  ii2-i  d£  (E-2) 

where  the  negative  sign  on  the  right  side  of  equation  (E-2) 
results  from  traveling  in  the  negative  -£  direction  when 
proceeding  from  point  1  to  point  3.   Assuming  a  very  small 
characteristic  mesh,  the  integrands  on  the  right  sides  of 
equations  (E-l)  and  (E-2)  can  be  taken  to  be  constant. 
With  this  assumption  and  some  rearranging,  the  equations 
become 

v3  =  7   (vj+vjj)  +  2-(e1-e2) 

l  sin  6  sin  Q 

+  j    [-sin  ^  __  A5l3  +  sin  ^  __£  Ar^] 

■*■  2 
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RIGHT- RUNNING 
CHARACTERISTIC 


STREAMLINE 


FIGURE  (E-l)   MODIFIED  CHARACTERISTIC  ELEMENT 


121 


and 


e3  =  i  ,v  -v2)  +  |  (61+e2) 


,  sin  6  sin  6 

+  7   [-sin  Ul  __  a^13  -  sin  y2  -— - 


An23] 
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Table  E-l   Results  of  Axisymmetric  M.  0.  C. 


POINT 

MACH  # 

y 

e 

V 

6+y 

e-y 

1 

1.082 

67.57 

91 

1 

23.4 

2 

1.256 

52.74 

95 

5 

42.3 

4 

1.435 

44.17 

100 

10 

55.83 

6 

1.605 

38.55 

105 

15 

66.45 

8 

1.775 

34.29 

110 

20 

75.71 

10 

1.95 

30.85 

115 

25 

84.153 

12 

2.134 

27.945 

120 

30 

92.05 

14 

2.14 

27.8 

122.65 

32.65 

94.85 

3 

1.198 

56.6 

90 

3.55 

146.6 

5 

1.35 

47.8 

93.57 

7.59 

141.37 

45.77 

7 

1.53 

41.0 

48.36 

12.62 

139.36 

57.36 

9 

1.70 

36.1 

103.23 

17.68 

139.33 

67.13 

11 

1.873 

32.31 

108.12 

22.72 

140.43 

75.81 

13 

2.05 

29.2 

113.03 

27.77 

142.23 

83.83 

15 

2.24 

26.5 

117.95 

32.83 

144.45 

91.45 

17 

2.349 

25.2 

120.58 

35.51 

145.78 

95.38 

16  1.49  42.1 
18  1.66  36.9 
20      1.835    33.0 


90.0 

11.65 

132.1 

47.9 

94.75 

16.71 

131.65 

57.85 

99.58 

21.81 

132.58 

66.58 
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Table  E-l   Results  of  Axisymmetric  M.  0.  C.  (Continued) 


POINT 

MACH  # 

y 

e 

V 

6  +  y 

e-y 

22 

2,02 

29.7 

114.43 

26.89 

134.13 

74.73 

24 

2.21 

26.9 

109.30 

31.98 

136.20 

82.40 

26 

2.415 

24.5 

114.18 

37.09 

138.68 

89.68 

28 

2.52 

23.3 

116.78 

39.80 

140.08 

93.48 

21 

1.835 

33.0 

90 

21.81 

123.0 

23 

2.02 

29.65 

94.82 

26.95 

124.47 

65.17 

25 

2.215 

26.85 

99.65 

32.08 

126.50 

72.80 

27 

2.42 

24.40 

104.50 

37.22 

128.90 

80.1 

29 

2.64 

22.20 

109.36 

42.40 

131.56 

87.16 

31 

2.77 

21.15 

111.95 

45.15 

133.10 

90.80 

34 

2.218 

26.75 

90 

32.18 

116.75 

36 

2.42 

24.35 

94.81 

37.36 

119.16 

70.46 

38 

2.65 

22.18 

99.64 

42.58 

121.82 

77.46 

40 

2.90 

20.10 

104.48 

47.86 

124.58 

84.38 

42 

3.045 

19.15 

107.06 

50.67 

126.21 

87.91 

39      2.655    22.13     90.0      42.63      112.13 

41      2.9      20.1      94.81     47.95      114.91    74.71 


124 


Table  E-l   Results  of  Axi symmetric  M.  0.  C.  (Continued) 


POINT 

MACH  # 

y 

e 

V 

e+y 

e-y 

43 

3.19 

18.25 

99.63 

53.37 

117.88 

81.38 

45 

3.36 

17.31 

102.20 

56.26 

119.51 

84.89 

18.23 

90.00 

53.42 

108.23 

16.45 

94.82 

59.03 

111.27 

78.37 

15.55 

97.38 

62.03 

113.93 

81.83 

52  3.19 
54  3.53 
56      3.73 


57      3.94     14.70     90.00     64.98      104.7 

59      4.19     13.82     92.55     68.16      106.37    78.73 

70      4.47     12.93     90.00     71.49      102.93 


72 

4.47 

12.93 

90.00 

71.49 

102.93 

77.07 

61 

4.22 

13.70 

92.10 

68.66 

105.80 

78.40 

58 

3.78 

15.35 

96.97 

62.68 

112.32 

81.62 

47 

3.41 

16.92 

101.78 

57.01 

118.70 

84.86 

44 

3.08 

18.95 

106.61 

51.41 

125.56 

87.66 

33 

2.8 

20.9 

111.47 

45.87 

132.37 

90.57 

30 

2.56 

23.0 

116.24 

40.49 

139.24 

93.24 

19 

119.92 

36.24 
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Table  E-l   Results  of  Ax i symmetric  M.  0.  C.  (Continued) 


POINT 

MACH  # 

y 

9 

V 

6+y 

e-y 

74 

4.47 

12.94 

90.00 

71.49 

102.93 

77.07 

63 

4.24 

13.65 

91.87 

68.87 

105.52 

78.22 

60 

3.79 

15.32 

96.72 

62.88 

112.04 

81.40 

49 

3.42 

17.03 

101.5 

57.18 

118.53 

84.47 

46 

3.1 

18.85 

106.31 

51.55 

125.16 

87.46 

35 

2.81 

20.8 

111.15 

45.97 

131.95 

90.35 

32 

115.90 

40.55 

76 

4.47 

12.93 

90.00 

71.49 

102.93 

77.07 

65 

4.27 

13.54 

91.51 

69.20 

105.05 

77.97 

62 

3.81 

15.22 

96.33 

63.15 

111.55 

81.11 

51 

3.43 

16.96 

101.08 

57.41 

118.04 

84.12 

48 

3.105 

18.79 

105.86 

51.73 

124.65 

87.07 

37 

2.82 

20.8 

110.67 

46.12 

131.47 

89.87 

78 

4.47 

12.93 

90.00 

71.49 

67 

4.315 

13.41 

91.03 

69.65 

104.44 

77.62 

64 

3.84 

15.11 

95.80 

63.59 

110.91 

80.69 

53 

3.455 

16.8 

100.49 

57.82 

117.29 

83.69 

50 

3.125 

18.63 

105.21 

52.10 

123.84 

86.58 

80* 

3.5 

99 

* 

estimated 
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APPENDIX  F 

A.   EJECTOR  NOZZLE  DESIGN 

To  assure  starting  of  the  cylindrical  nozzles  a  back 
pressure  must  exist  which  is  lower  than  the  total  pressure 
downstream  of  the  shock.   For  a  stagnation  pressure  of  ten 
atmospheres,  the  total  pressure  downstream  of  a  normal 
shock  at  a  Mach  number  of  five  is  less  than  one  atmosphere. 
Ejectors  were  designed  to  lower  the  back  pressure  to  a  value 
below  the  total  pressure  downstream  of  the  shock  assuming 
Mach-five  flow. 

The  models  used  for  finding  the  flow  conditions  at  the 
ejector  nozzles  are  shown  in  figures  (F-l)  and  (F-2) . 
Applying  the  continuity  equation 

PlUlAl  +  p2U2A2  =  P3U3A3   '  (F-1) 

the  momentum  equation 


P1U12A1  +  P^  +  P2U22A  +  P2A2  =  P3U32A3  +  P3A3  ,  (F-2) 


and  the.  energy  equation 


P1U1A1(CpT1  +U12/2)  +  p2U2A2(CT2  +U22/2)  =  P^^  (CT3  +  U"32/2) 

(F-3) 
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FIGURE  (F-l)   EJECTOR  NOZZLE  MODEL 
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FIGURE  (F-2)   TOP  VIEW  OF  EJECTOR  FLOW  MODEL 
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to  the  control  volume  shown  will  yield  the  flow  conditions 
necessary  at  station  one  to  assure  the  desired  back 
pressure. 

The  following  assumptions  are  made: 

(1)  isentropic  flow  from  station  4  to  2 

(2)  isentropic  flow  from  reservoir  to  ejector  nozzle  exit 

(3)  P3  =  1  atm. 

(4)  PT   =  9.5  atm. 

(5)  P2  -  Px 

(6)  TTl  =  520°R 

(7)  M  =  5.0  at  shock 

Using  normal  shock  relations  the  known  flow  at  station  (5) 
can  be  used  to  obtain  conditions  at  station  (4) .   For 
M  =  5.0,  the  pressure 


P4/P5  =  29.0 

P4  =  (29)  (39.99) 
P.  =  .548  atm. 


density 

P4/P5  =  5.0 

p4  =  (5.0) (pQ) (p5/pc) 

p4  =  (5.0)  (.02378  g)  (.01134) 

p4  =  .0434  lbm/ft3  , 
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and  temperature 


T4/T5  =5.8 

T4  =  (5.8) (TQ) (T/TQ) 
=  (518)  (520)  (.1667) 


T4  =  502.76  °R 


can  easily  be  found. 

Using  the  isentropic  relations 


P4/P2  =  (p2/P4)  =  (T2/T4)   ' 


all  conditions  at  station  two  can  be  found  in  terms  of  the 
static  pressure: 


and 


p2  =  p4(Vp2)1/y 


p2  =  (.0434)  (J^i8-)1/1-4 


T2  =  T4(VP2}  (Y"1)/Y 


T2  =  (502.76) ( .548)4/1.4 

P2 


•     • 


Since  m2  =  m4  ,  the  velocity  at  2  can  also  be  found 


m 


2  =  P2U2(2A2)  =  A4 
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where  the  area  is  taken  as  (2«A2)  since  the  control  volume 
inlet  area  (A~)  is  only  one-half  of  the  total  exit  area  of 
the  flow.   Rearranging: 


m, 


ut  =  htta k — rr  *  (144) 

2    p2  [2  (A3  -  A1)  ] 


(6.463)  (144) 


p2(2)  (9  -  .785) 


U  =  55'9 


2     P2 


All  conditions  at  station  two  are  now  known  in  terms  of  P2 . 
Again  using  isentropic  relations,  conditions  at  station  one 
can  be  determined  in  terms  of  P, : 


PTn  -(Y-D/Y 
Tl  =  TT  (TJl) 


o  5  -.4/1.4 
T,  =  520  (^2.)  °R 

1 

and  PT 

pl  =  p0  ("pf  > 


9  5  -1/;U4         3 
P;L  =  (.02378  g)  (^2.)         lbm/ft 
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To  find  the  velocity  at  station  one 


T   =  T  (1  +  *^,   M  ) 
iT1     1^  ^T^     1  ; 

2    "l2 
Since   M, 


1    yRT,   ' 


2yRT, 
Ul  =  V  ^  ~  1}  (TrTL)  g       ft/sec 


Hence  all  conditions  at  station  one  are  known  in  terms  of 

pi- 

If  P,  (=P2)  ^s  known,  equations  (E-l)  and  (E-2)  can  be 
combined  to  find  IU  and  p.,.   A  computer  program  was  written 
to  vary  the  values  of  D-,  and  D^  and  to  iteratively  increase 
P,  from  a  set  value  until  it  satisfied  the  energy  equation. 
An  explanation  of  the  program  EJECTOR,  a  listing,  and  the 
results  in  tabular  form  are  found  in  this  appendix.   Figure 
(26)  shows  the  results  in  graphical  form.   The  broken  line 
at  a  pressure  of  0.617  represents  the  total  pressure  down- 
stream of  the  shock  assuming  P   (upstream)  =  10  atm.  and 
M  =  5.0.   Acceptable  combinations  of  D,  and  D3  must  result 
in  a  pressure  P~  below  this  line. 
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B.   EXHAUST  LENGTH 

The  length  of  the  exhaust  stack  must  be  long  enough  to 
ensure  complete  mixing  of  the  primary  flow  from  the  ejector 
nozzles  and  the  secondary  flow  from  the  cylindrical  nozzles. 

The  model  used  for  mixing  of  the  two  flows  is  shown 
in  figure  (F-3) . 

Bailey  [Ref.  6]  defines  a  non-dimensional  spread 
parameter  which  is  a  good  approximation  to  experimental  data 


a  =  24  ^H  (F-4) 

ul 


where 


Ul+U2 
U      =     ~ 

mean      2 


The  non-dimensional  similarity  variable  for  the  mixing 
zone  was  defined  to  be 


n  =  —   .  (F-5) 


A  value  of  r\   =  1.84  gives  the  outer  edge  of  the  mixing  zone 
where  the  velocity  equals  99%  of  the  centerline  value. 

The  minimum  value  of  the  mixing  length  is  where  the  edge 
of  the  mixing  zone  touches  the  side  walls,  or  when 


D-D 

Y  =  Y     =    * 
max      2 
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FIGURE  (F-3)   EJECTOR  FLOW  MIXING  MODEL 
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Substituting  Y  =  Y    and  equation  (F-4)  into  equation  (F-5) , 
J  max  .  ' 

the  minimum  mixing  length.  X  .  ,  can  be  found: 

3     =»   *   mm 

_        max 
min     r\ 

U       (D,  -  D,  ) 

24   mean   v  3    1' 

~U^    2(1.84) 


(U,  +u9) 
6.522(4.85-1.75) 


2U1 


"min'^  (1  +  07)  (3'1J    ■  (F"6) 


Using  values  from  program  EJECTOR  for  the  combination  of 
D-,  and  D3  chosen: 


U1  =  1860  ft/sec 
U2  =   664  ft/sec   . 

Putting  these  values  in  equation  (F-6) : 

Xmin  =  i4^(1+&  13.1) 

X  .   =  13.72  inches 
min 

The  exhaust  stack  must  therefore  extend  a  minimum  of  13.72 
inches  beyond  the  ejector  nozzle  exit.   Figure  (F-4)  shows 
the  length  of  the  exhaust  stack  with  the  ejector  nozzles 
installed. 
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C.   NOZZLE  SHAPE 

For  the  chosen  combination  of  D,  =  1.75  inches  and 
Do  =  6.5  inches 


U,  =  1860  ft/sec 


(See  section  D  of  this  Appendix.)   This  velocity  is  clearly 
supersonic.   To  determine  the  throat  area  required  for  the 
given  exit  diameter  and  flow  velocity,  the  exit  Mach  number 
must  first  be  found.   Since 


2 

M   =  _ 

A^ 


m2  -  U2 

M   = 


yRT    ' 


and 


T  =  T/(l  +^5^.  M2)   , 


then,  after  substituting  and  rearranging,  the  Mach  number 
can  be  expressed  as 


VI  -  IzL  (K) 


where 


U2 
K  = 


YRgTQ 
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Inserting  the  appropriate  values 


m  -  -a/  2.77 
M  "  VI -.554 


M  =  2.49 


For  isentropic  flow  through  the  nozzle 


Y+l 
A/A*  =  J  [(  ^  )(i+I-iM2)]  TFFTT  ^ 


M  L  VY+T 


For  M  =  2.49 


A/A*  =2.61   . 


Solving  for  throat  height: 


A A (2T6T} 


A*  =  .921   in2   , 


or 


/4A* 

D*  =  Vir 


D*  =  1.08  inches 
The  ejector  nozzles  were  fabricated  from  plexiglas  with  the 
throat  diameter  of  1.08  inches.   Figure  (F-4)  shows  the 
details  of  the  ejector  nozzle  geometry. 
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FIGURE  (F-4)   EJECTOR  NOZZLE  GEOMETRY 
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D.   EJECTOR 

The  program  EJECTOR  solves  the  continuity  and  momentum 
equations  for  different  combinations  of  nozzle  diameter  and 
exhaust  stack  width  assuming  a  starting  value  of  pressure 
at  stations  one  and  two  of  0.37  atmospheres.   The  pressure 
P,  is  increased  in  increments  of  0.0001  atmospheres  until 
the  flow  conditions  at  1,  2  and  3  satisfy  the  energy  equation. 
When  the  energy  equation  is  satisfied,  the  values  are  printed, 
and  the  process  is  repeated  with  new  values  of  D-,  and  D,. 

Inputs  to  the  program  include  conditions  immediately 

3 

downstream  of  the  shock:   p.  (lb  /ft  )  ,   T,  (°R)  , 

P.  (atm.)   and  m,  (lb  /ft  )  ;   and  total  temperature  (°R) 
4  4    m 

and  total  pressure  (atm.)  at  station  one. 

Outputs  are  D-,  and  D^  (both  in  inches)  ,  P,  (=P2)  in 
atmospheres,  and  U,  and  U2  ,  both  in  ft/sec. 
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Table  F-l   Symbols  for  Program  EJECTOR 


Computer  Notation   Algebraic    Description 

2 

Flow  area  (in  ) 

Static  pressure  (atm.) 

Static  density  (lb  /ft3) 

Mass  flow  (lb  /sec) 

Specific  heat  ratio 

Gas  constant  for  air 
(53.3  ft-lbf/lb   °R) 

Specific  heat  @  constant 
pressure 

778.26  ft-lbf/BTU 

Static  temperatures  (°R) 

Stagnation  temperature  (°R) 

Stagnation  pressure  (atm.) 

Velocities  at  stations  1,  2 
and  3  respectively  (ft/sec) 


A  A2  A3 

Al 

A2 

A3 

PI  P2  P3 

P4 

Pl 

P2 

P3 

P4 

Rl  R2  R3 

R4 

Pl 

P2 

p3 

p4 

Ml  M2  M3 

M4 

mx 

m2 

m3 

&4 

G 

Y 

R 

R 

CI 

cp 

J 

j 

Tl  T2  T3 

T4 

Ti 

T2 

T3 

T4 

T 

TTl 

P 

PTl 

Ul  U2  U3 

Ul 

U2 

U3 

K 

2*g 

p 

L 

C  *J 
P 

Dl 

D, 

Diameter  of  one  of  the  two 
ejector  nozzles  (inches) 

D3  D^  Length  of  one  side  of  exhaust 

stack. 


141 


D.       LIST    OF   COMPUTER   PROGRAM    "EJECTOR"    AND    RESULTS 


'R1N1      '  IHPUl       '  il  !HIT  EONS  S  ! 

,    j         '      |    |-,]|     j      ■  •      j|J     .       J     : 

M-?l      .M4 
.  ■  i ;         N  iTRL        i {■  I D I      0! 

P !     N'i     "  J.  H  r  U  i     Ml  i  R )  ?  P  I  1  k  H  i  I'l  / 

;    ■  i 

I  48    J"™ 

1  58    L  '  ■■  0 

1  68    PR  IN        [i !  "  i     D3 ::  ?  "PI !'    ::  U i "  s  "  U2 

I     I'''  !'J         I      :  '      [Hi 

1 ??    P1=0. 37 

1  f-'  o     r  U  K     J.J  '.:.' ::":  fa      i  U     3     o  i  E  P     U  „  b 

180    R3=D3  i\2  ••••'!  44 

;  36    R2=fi3  -fl'   (  3*  (.  D3  '3  )  £B  J  ••''3.  44 3--  2  ) 

219    R2=R4  '  (P4/P2)  N  1/G) 

£30  U2=M4/<R2*fl 

-    .  1 1  i     K  1 :::: '-.  P  ■■ '' '•  R " v  !  )  ■'  '■  '•■  P  ■■' '  P  1  ,';  '*'  '•  ' "'  I  ■"''  I .:;  -;  *  2  1 

■/  f. j'-'i       j  !  1    ::: ;  g  Q  p  (   (   (    J  ., '  "j"  j    )  j    )  *  V  #  fj  ¥•  P  ft '  |"   j   ■¥•  3  £ 

i:!  ■•''  0  i. ' ..:' ":::  k  1 *  U  1  *  H  •••''  H  3  +  K ,. :!  *  I  J  i '  #  1 1 2  -  •"'  hi  3 

230  ij3=  (i?  ]  ¥.|j  1 12  :tR/R3  *• '  P2  -F  3  )  ;:  21 1  fa 

2  9  0  R '' : :::  R  i  *  U  1  *  fl  •■•"'  (  U  3  *  R  3  ) •+  R  2  *  U  2  *  R  2  ■••'  ( 
333  T3=P3*21 16- <R*R3> 

313    N3 =R1*U1*R 
333    N2=R2*U2*R2 
333    M3=R3*U3*R  ■ 

i    L=C1 
333    Wi=Ml*<L*T3  ■    :  i      '    K) 

3  r'  "J     N ,: .:!  -::  M .:..:  '*  v  L  *  i  2  +  U  2 ''''' ,::!  /  K  ) 

333  N=Wi +W2 

:  •  !=M3*<1      !     HJ3t2    K) 

4Q@  D=kl->! 

::  IF  D<10  THEN  410 

J  L'l     '        T  C        :":   |        :  H       ..'.'  ■■■■'-        r  :  !  j.."  M        ■■   "   :  - 
"r  t.14       I"   1  —  I'"    !.  +  hj  „   r.M:;i  | 

A  0  fa    G  0  T  0    2  8  0 

4  i  o    PRINT    IU  j-D3»  PI  p  Ul  ?  U2 
4  2  0    r !  E  :A  T    D  3 

430  PRINT 
448  NEXT  HI 
458  END 
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